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ABSTRACT 

(Distribution  Limitation  Statement  No.  2) 

A  comparison  between  tvo  one-dimensional  Lagrangian  hydrocodes  has  been  made. 

The  two  hydrocodes  are  a  von  Neumann-Rich tmyer  hydrocode  (AFWL's  PUFF)  and  a 
Lax-Wendroff  hydrocode  (the  two-step  version  with  artificial  viscosity).  The 
comparison  was  made  by  applying  the  hydrocode  test  problems  as  described  in 
HYDROCODE  TEST  PROBLEMS,  AFWL-TR-67-127 ,  February  1968.  The  most  apparent 
difference  between  the  von  Neumann-Rich tmyer  hydrocode  and  the  Lax-Wcndroff  is 
the  greater  tendency  of  the  Lax-Wendroff  scheme  to  oscillate.  In  those  flows 
in  which  there  are  no  strong  shocks  or  strong  rarefactions  or  vacuums,  the  Lax- 
Wendroff  scheme  is  more  accurate.  However,  in  those  flows  in  which  there  are 
strong  shocks  or  strong  rarefactions  or  vacuums  the  von  Neumann-Rich tmyer  scheme 
is  more  accurate.  The  Lax-Wendroff  scheme  cannot  handle  vacuums  because  of  the 
use  of  the  specific  volume  instead  of  the  density  as  a  fluid  variable.  It 
appears  that  it  might  be  possible  to  combine  the  better  features  of  the 
von  Neumann-Richtmyer  and  the  Lax-Wendroff  schemes  to  produce  a  better  hydrocode. 
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SECTION  I 
INTRODUCTION 


This  report  is  the  description  of  a  comparison  between  AFWL's  PUFF  hydrocode 
and  LAX-WENDROFF  two-step  method  with  viscosity.  The  basis  for  comparison  is 
the  solutions  these  hydrocodes  produce  to  a  series  of  hydrocode  test  problems. 

The  problems  involve  shocks  and  rarefactions  and  interactions.  The  hydrocode 
test  problem  solutions  are  Known  exactly.  Brief  descriptions  of  these  test 
problems  will  be  given  here.  For  more  details  see  "Hydrocode  Test  Problems" 
AFWL-TR-67-127. 

1.  THE  PUFF  HYDROCODE 

Let  the  points  of  a  rectangular  network  with  spacings  Ax  and  At  be  denoted 
by  *£»  t  »  (Z.  m  03  1,  2,  •  ••.  Ls  n  *  0,  1,  2,  •••).  There  will  also  be  occasion  to 
deal  with  intermediate  points,  having  coordinates  x  ,  =  h  (y  +  x  ) ,  tn+?s  - 

,  1  £*r'$  £+1  £ 

*s(tn  +  tn).  To  facilitate  the  writing,  introduce  abbreviations  such  as 
v”  +■, 5  c“>-  *«• 


Un-M  _  yn— 

Z _ £ _ 

At 


z+h 


n-h  n  n-*s 


(ZW 


mz-H)/2 


(1) 


is  PUFF's  difference  approximation  to  pQ  |£  -  -  where  ZM  is  the  zone  mass, 

U  is  the  fluid  velocity,  P  is  the  fluid  pressure  and  q  is  the  artificial  viscosity. 


Yn+1  yn 

h _ h 

At 


,n+^ 


(2) 


is  PUFF's  difference  approximation  to  =  U,  where  X  is  fluid  position. 


I 
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zm„  , 

Dn+1  - - - 

l-h  j.n+1  -  xn^ 
Z  ~  l- 1 


p0  3X 


is  PUFF's  difference  approximation  to  ,  where  p  is  the  fluid  density.  Now 

let  4U  -  uf*  -  U £*  . 


Then  PUFF's  q  is  given  by 


■co  -  C1 


/.n+]  ,  n  \ 

.csn-l)lu 


where 


C0  -  1.8 


Cj  -  .25 


CS  -  isothermal  sound  speed 


CS2  if 

dp 


e  const. 


e  is  specific  internal  energy. 


en+l  .  „n  fpn+l  +  n+4  +  pn  +  qr.-js\ 

o  m  +  Vi-H  qi-h  Jdd  . 


is  PUFF's  difference  approximation  to 


0  -  +  (P+q)  ^  f 


l 


2 


*3 
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which  results  from 


n  3e  .  ,D.  ,3V  ,3V  1  30 

0  =  It  +  (p+q)  It"  and  I?  "  lx" 


where  V  *  —  is  the  specific  volume. 


Lastly,  the  equation  of  state: 


pn+l  p  /  n+1  n+l\ 

P£-*1  Pl*W  °a-h) 


(6) 


PUFF's  method  of  solution  is  this:  Suppose  all  quantities  are  known  for 

superscript  n  or  n-h  (this  is  referred  to  as  being  at  cycle  n) .  Compute 

for  each  £  from  (1),  then  compute  for  each  £  from  (2),  then  compute  pj^ 

for  each  £  from  (3),  then  compute  q^+?  for  each  £  from  (A),  then  compute  e”+£ 

for  each  £  by  simultaneously  solving  (13)  and  (14).  At  this  point  ail  variables 
have  been  advanced  to  cycle  n+1.  Next  PUFF  does  its  time-step  computation. 


At  *  .9  min 


Yn+1  Yn+1 

X£  h-1 


£  CSJ^(1+2-Cj)-.4C2aU 


(7) 


where,  as  in  the  q  calculation. 


1.8,  C,  =  .25 


Remember 


CS‘ 


dP 

dp 


e  const. 


Therefore  CS  is  the  isothermal  sound  speed. 
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The  isentropic  sound  speed  is  defined  by 


C2 


S  const. 


where  S  is  the  entropy. 

For  a  y  ~  law  gas  C2  =  yCS2.  Therefore  for  AU  very  small 


At  « 


Q  ,  X«+l  -  x?+* 
.9  y  min  _J. _ £-1 


(3/2)  * 


-n+1 


7572)  * •?  for  y  * 


if 


At  «  e 


yn+l  _  yn+1 
rain  l  l-l 

l  rn+l 


0  is  called  the  effective  CFL  number. 

For  further  details  about  PUFF  see  AFWL-TR-66-48  and  AFWL-TR-67-127. 
2.  THE  LAX-WENDROFF  METHOD 


The  LAX-WENDROFF  two-step  method  with  viscosity  uses 


■  ’>(“"« +  °")  -  (ft + ")(f;+1  -  f") 


and 


Un+1  =  un  -  —  fFn+|i  - 

.1  j  4z  \  j+4  i-h) 


(8) 


(9) 


as  the  difference  approximation  to 


3U  5F(U) 

3t  dz  ' 


4 
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q 


(10) 


where  d  =  Vqc/V,  c  is  isentropic  sound  speed,  V  is  specific  volume  and  VQ  is  a 
constant  with  dimensions  of  specific  volume  defined  by 


*  =  v0  yp0(x)dx 

In  our  case  VQ  =  1  therefore,  z  is  the  Lagrangian  mass  variable,  b  is  a  dimen¬ 
sionless  parameter  which  was  chosen  to  be  .5. 


U  = 


where  E  *  e  +  *5u2  and  e  is  the  fluid's  specific  internal  energy,  u  is  the  fluid 
velocity,  and  P  is  the  fluid  pressure. 

The  time  step  restriction  is 


At 


4*rM) 


min 


xi  -  xm 

V* 


where  X  is  fluid  position. 


For 


b,%,(1  +  £f-  t. . 


78 


LAX-WENDP.OFF ' s  method  of  solution: 

Suppose  all  quantities  are  known  for  superscript  n.  Compute  for  each  j  from 

(8)  then  compute  for  each  j  from  (9).  Now  all  variables  are  advanced  to 

cycle  n+1.  Next  LAX-WENDROFF  does  its  time  step  computation 


5 


I 


If 

% 

f 

i 

I 

i 

! 

i 

k 


l 
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where  0  -  .78. 


At 


0 


min 

2. 


Yn+1 


Yn+1 

X£-l 


rn+l 

C£-Js 


For  more  details  see  Richtmyer  and  Morton:  Difference  Methods  for  Initial  Value 
Problems,  Interscience  Publishers,  a  division  of  John  Wiley  and  Sons,  1967. 

3.  THE  HYDROCODE  TEST  PROBLEMS 

Since  the  geometry  is  one-dimensional  slab,  the  problems  may  all  be  thought 
of  as  flows  in  a  smooth  pipe  of  constant  cross  section.  There  are  seven  problems. 
The  first  problem  is  the  flow  that  results  from  a  piston  moving  into  the  gas  with 
a  constant  velocity.  The  second  problem  is  the  flow  that  results  from  pulling  a 
piston  away  from  the  gas  with  a  constant  velocity.  The  third  problem  is  the  flow 
that  results  from  a  piston  moving  into  the  gas  with  a  constant  acceleration.  The 
fourth  problem  is  the  flow  that  results  from  pulling  the  piston  away  from  the  gas 
with  a  constant  acceleration.  The  fifth  problem  is  the  flow  that  results  by 
removing  a  partition  between  two  different  states  of  the  gas  at  rest.  The  sixth 
problem  is  the  flow  that  results  from  the  collision  of  two  shock  waves.  The 
seventh  problem  is  the  flow  that  results  when  one  shock  overtakes  another  one. 

For  more  details  see  AFWL-TR-67-127. 
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SECTION  II 

COMPARISON  OF  THE  HYDROCODES 

For  each  test  problem,  the  exact  solution,  the  PUFF  solution,  and  the  LAX- 
WENDROFF  solution  will  be  described.  In  describing  the  PUFF  and  LAX-WENDROFF 
solutions  an  error  tabic  and  graphs  of  their  solutions  will  be  used  to  compare 
with  graphs  of  the  exact  solutions.  In  the  error  table  the  numbers  presented  are 
labeled  Sum  Abs.  Error,  Sum  Sqr.  Error,  and  Maximum  Error.  These  numbers  are  now 
defined.  Let  Pp(J)  be  the  PUFF  pressure  in  zone  J  and  let  PE(J)  be  the  exact 
pressure  in  zone  J.  Let  PM  be  the  maximum  of  the  P£(J). 


Sum  Abs.  Error  (for  P) 


I 

J 


PP(J)  -  PE(J)| 


rM 


Sum  Sqr.  Error  (for  P) 


Maximum  Error  (for  P) 


mnx|Pp(J)  -  PE(J)| 

- - ^ - sgh(pp<jm)  -  PEUM)) 


JM  is  the  zone  index  of  the  maximum  error  and  SGN  is  the  sign  function.  The 
error  functions  are  likewise  defined  for  the  velocity,  density,  and  energy 
(specific  internal). 

Also  in  the  error  table  are  presented  the  sums  of  the  internal  energy,  kinetic 
energy,  and  total  energy  of  the  exact  solution,  FUFF  solution,  and  LAX-WENDROFF 
solution  (the  unit  is  ergs).  In  addition,  the  error  table  contains  the  problem 
time,  computer  time  (CP  time  on  the  CDC  6600),  cycle  number,  and  the  number  of 
active  zones. 

The  graphs  are  organized  in  this  manner:  pressure  and  density  are  plotted 
in  the  same  graph  as  are  velocity  and  energy  (specific  internal). 


7 
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1.  TEST  PROBLEM  SCTP-I 
a.  The  Exact  Solution 

In  this  problem  a  piston  moves  to  the  right  into  the  gas  at  a  constant 
velocity.  The  solution  has  a  steady  profile.  The  solution  profile  is  two  con¬ 
stant  states  separated  by  the  shock  discontinuity.  That  is,  each  fluid  parameter 
(pressure,  density,  fluid  velocity,  etc.)  is  a  constant  from  the  piston  face  to 
the  shock  and  each  is  another  constant  to  the  right  of  the  shock.  The  symbols 
used  to  describe  the  problem  further  are 

sound  speed  to  the  left  of  the  shock 

Cr  sound  speed  to  the  right  of  the  shock 

P^  pressure  to  the  left  of  the  shock 

Pr  pressure  to  the  right  of  the  shock 

density  to  the  left  of  the  shock 

Pr  density  to  the  right  of  the  shock 

specific  volume  to  the  left  of  the  shock 

Vr  specific  volume  to  the  right  of  the  shock 

v^  fluid  velocity  to  the  left  of  the  shock 

Vp  piston  velocity 

vf  fluid  velocity  to  the  right  of  the  shock 
Vg  shock  velocity 
piston  position 


~P 

*Q 


quiet  zone,  a  position  far  enough  to  the  right  so  that  the  f.-s  is 
still  at  rest;  energy  sums  are  taken  out  to 


Xg  shock  position 

There  are  two  variations  of  SCTP-i  and  these  are  denoted  SCTP-I-A  and  SCTP-I-B. 
For  SCTP-I-A  the  piston  is  started  at  C,  and  the  shock  is  started  at  50  meters, 
with  the  fluid  parameters  on  the  right  at 

Pr  *  104  dynes /cm2 

p^  *  10-6  gm/cm^ 

vr  «  0.  cm/sec 


I 


MM 
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Xq  *  300  meters 
This  yields 

C  *  /y  x  105  cm/sec,  which  for  y- 1.4  yields 
Cr  =  1.18  x  1G5  cm/sec 
Then  one  sets 

Vp  ■  C  =  1.18  x  10s  cm/sec  and  this  yields 
v  =  1.18  x  105  cm/sec 
vg  =  2.09  x  10®  cm/sec 
P  -  3.47  x  104  dynes /cm2 
V,  =  4.34  r.  10®  cm3/gn 

This  problem  is  run  for  .1  second,  with  initial  zones  of  1  meter.  For  SCTP- 
I-B  the  shock  is  again  started  at  50  meters  and  with  the  fluid  parameters  on  the 
right  at 

Pr  =  104  dynes /cn2 
Pr  =  10” 6  gm/cm3 
vr  =  0.  cm/sec 
Xq  =  300  meters 

This  yields 

Cr  =  1.18  x  10s  ca/sec 
Then  one  sets 

Vp  =*  100  Cr  »  1.18  x  107  cm/sec  and  this  yields 
v.  =  1.18  x  107  cm/sec 

L 

P,  -  1.68  x  10®  dvnes/cm2 
V  =  1.67  x  10®  cm3/ga 
C.  =  6.26  x  10®  cm/sec 

L 

vs  =  1.42  x  10'  cm/sec 

This  problem  is  run  for  L0~3  seconds  with  initial  zones  of  1  meter. 
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b.  The  PUFF  Solution 

The  regions  where  the  largest  errors  occurred  were  the  regions  where  the 

shock  was  initially  and  where  the  shock  is  currently.  See  Tables  I-A  and  1-B 

and  Figures  I-A  and  I-B. 

c.  The  LAX-WENDROFF  Solution 

SCTP-I-A  was  run  with  a  viscosity  facto::  of  .5  and  a  time  factor  of  .78. 
aCTP-I-B  was  run  with  a  viscosity  factor  of  .5  and  a  time  factor  of  .25.  In 
order  to  get  SCTP-I-B  to  run  it  was  necessary  to  start  off  with  10  time  steps 
with  ^ero  viscosity  factor  and  .025  time  factor. 

The  regions  where  the  largest  errors  occurred  were  the  regions  where  the 

shock  was  initially  and  where  the  shock  is  currently.  The  main  difference  between 

the  PUFF  and  LAX-WENDROFF  solutions  is  the  oscillations  behind  the  shock  front. 

The  oscillations  are  much  more  pronounced  in  the  LAX-WENDROFF  code.  See  Figures 
I-A  and  I-B  and  Tables  I-A  and  I-B. 
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Sum  Int.  Energy  Sum  Kin.  Energy  Sum  Tot.  Energy 

EXACT  3.09498  x  1012 _ 3.09324  x  1012 _ 6.18823  x  1012 

LAXWEN  3.09016  x  1012  3.08281  x  1012  6.17297  x  1012 


EXACT  SOLUTION  TIME  1-000*10 


RADIUS >10  (CM) 


EXACT  SOLUTION  TIME  l.Q00*10 


RAOIUS *10  (CM) 


•Till 


PO-PUFF 


CYCLE  1033 
PUFF  66  mm  SCT 


RADIUS  *10  ’(CM) 
Figure  I -A.  VE-PUFF 


SMS 


igure  I -A, 


CYCLE  1463  TIME  1.000*10 

PUFF  66  »»  SCTP-I-B 


TIME  1 .000*10 
IflXWEN 


•LAX-WENDROFF 


riME  1.000*10 
LflXWEN 


|-B.  VE-IAX-WENDROFF 
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2.  TEST  PROBLEM  SCTP-II 
a.  The  Exact  Solution 

In  this  problem  a  piston  pulls  away  with  constant  velocity  from  the  gas 
at  rest  in  a  pipe.  The  piston  moves  to  the  left  with  constant  velocity  vp  <  0 
away  from  the  gas  on  the  right.  This  causes  a  rarefaction  wave  to  move  to  the 
right.  For  a  graphical  description  see  the  exact  solution  plots  in  Figures  II. 
The  exact  solution  for  the  velocity  is  piecewise  linear  as  a  function  of  X. 
Starting  at  the  piston  on  the  left  at  position  Xp(t)  the  velocity  is  the  constant 
Vp  from  Xp(t)  to  what  is  called  the  back  of  the  rarefaction  wave  and  denoted 
XR(t).  From  XR(t)  rightwards  to  X^(t)  the  velocity  rises  from  vR  linearly  to 
zero.  Xc(t)  is  the  front  of  the  rarefaction  wave.  To  the  right  of  X^(t)  the  gas 
is  at  rest  so  the  velocity  is  a  constant  zer?. 

Xp(t)  =«  Xp(0)  +  vpt 

XR(t)  «  Xp(0)  +  (cr  +  vp)t 

Xc(t)  -  Xp(0)  +  crt 

The  rest  of  the  variables  are  then  determined  by  the  simple  wave 
formulas: 


C(X, t)  =  cr  +  ^  v(X,t) 


There  are  five  variations  on  this  problem.  This  much  is  common  to  all 

of  them: 

Pr  =  ]04  dynes/cm2 
Pr  *  10*’1’  gm/em-5 
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Cr2  «  YPr/or  ■  1.4  x  1010  cm2/oec2 
AX  =  100  cm 
Xp(O)  =  100  meters 

Xq  “  300  meters 

and  all  variations  are  run  out  to  .1  second.  The  variations  are  in  the  piston 
velocity. 

SCTP-II-A  J vp |  -  Cr/(>+l) 

SCTP-II-B  jvp{  -  2Cr/(y+l) 

SCTP-II-C  jvp |  -  2Cr/(y-l) 

SCTP-II-D  |vp|  -  4Cr/(y-l) 

SCTP-II-E  Free  boundary  condition  on  the  left  in  place  of 

withdrawing  piston  condition.  That  is,  it  is  as 
if  at  time  zero  one  removes  a  separator  to  the 
left  of  which  is  a  vacuum. 

These  variations  were  introduced  to  investigate  the  codes  response  to 
the  following  situations:  in  A,  XR(t)  moves  to  the  right  with  velocity  Cr/2;  in 
B,  XR(t)  is  stationary.  In  both  A  and  B  the  piston  is  not  pulled  out  too  fast 
for  the  gas  to  follow;  therefore  the  pressure  and  density  are  positive  constants 
from  Xp(t)  to  XR(t).  However,  in  C,  D,  E,  XR(t)  moves  to  the  left  with  velocity 
-2Cr/(y-l)  and  between  Xp(t)  and  XR(t)  there  is  a  vacuum.  In  C  the  piston  is 
pulled  out  with  exactly  the  escape  velocity  of  the  gas,  -2Cr/(y-l),  therefore, 
XP(t)  »  XR(t).  In  D  the  piston  is  pulled  out  faster  than  the  gas  can  follow  and 
so  Xp(t)  <  XR(t) •  In  E  the  code  is  allowed  to  compute  its  own  escape  velocity. 

b.  The  PUFF  solution 

On  A  and  B  PUFF  tended  to  underround  at  Xq  then  overround  at  XR  and 
undershoot  just  to  the  left  of  XR.  See  Tables  and  Figures  II-A  and  II-B.  In 
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C,  D,  E  PUFF  again  tended  to  underround  at  X^.  In  E  the  gas  front  did  not  move 
as  far  to  the  left  as  it  should.  This  is  because  of  the  finite  mass  in  the  left 
hand  zone.  If  the  zoning  were  made  finer  to  the  left  so  that  the  left  hand  zone 
had  a  smaller  mass,  then  the  left  hand  zone  would  move  out  mure  nearly  at  the  rate 
at  which  the  gas  should  escape.  See  Tables  and  Figures  II-C,  -D,  ~E. 

c.  The  LAX-WENDROFF  Solution 

Because  the  LAX-WENDROFF  scheme  uses  specific  volume  instead  of  density 
it  is  able  only  to  run  SCTP-II-A  and  B.  This  is  because  of  the  vacuums  in  C,  D, 
and  E.  In  the  vacuum  the  density  is  zero  and  the  specific  volume  is  infinite. 

If  this  scheme  is  to  be  used  for  problems  in  which  there  are  vacuums  or  near 
vacuums  the  specific  volume  must  be  changed  over  to  density.  The  time  step 
factor  used  was  .78  and  the  artificial  viscosity  factor  used  was  .5. 

In  A  and  B  the  LAX-WENDROFF  scheme  tended  to  underround  at  Xq,  then  over¬ 
round  at  XR,  then  undershoot  to  the  left  of  XR,  and  then  dampingly  oscillate 
toward  the  left. 

Lastly,  one  zone  to  the  right  of  the  piston  face  there  is  a  density  dip. 

It  is  believed  that  the  reason  the  LAX-WENDROFF  scheme  has  larger  errors  on 
this  problem  is  because  it  has  a  q-factor  (artificial  viscosity)  in  expansion. 

In  the  PUFF  code  q  is  not  used  in  expansion — only  in  compression.  It  seems  that 
for  the  LAX-WENDROFF  scheme  to  compete  with  PUFF  its  q  needs  to  be  modified  also 
for  compression  as  noticed  in  SCTP-I.  One  other  thing  at  this  point:  notice 
that  the  computer  times  for  the  LAX-WENDROFF  scheme  are  longer.  This  is  because 
PUFF  is  a  production  code  and  much  time  was  spent  making  it  run  efficiently.  On 
the  other  hand  no  time  was  spent  trying  to  make  the  LAX-WENDROFF  scheme  coding 
efficient.  The  difference  scheme  was  merely  programmed  to  test  its  accuracy. 

So,  if  the  LAX-WENDROFF  scheme  is  to  be  used  for  production,  then  time  should  be 
spent  in  making  the  programming  more  efficient. 
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_ Sum  Int.  Energy  Sum  Kin.  Energy  Sum  Tot.  Energy 

EXACT _ 4.43248  x  108 _ 2.92345  x  107  4.72482  x  108 

LAXWEN  4.43411  x  108  2.89319  x  107  4.72343  x  108 


LAXWEN 


ERRORS  ON  SCTP-II- 


Int.  Energy  Sum  Kin*  Energy  Sum  Tot*  Energy 


ERRORS  Oil  SCTP-II-E 


LAKWEN 


EXPCT  SOLUTION  TIME  1 .000*10 
SCTP-1I-P _ 


PP01U5  UQ  (CM  I 
Figure  1 1- A.  PD -EXACT 


EXRCT  SOLUTION  TIME 
SCTP-II-R 


Figure  1 1 -A.  VE -Exact 


SCTP-I I-fl 


Figure  1 1 -A.  VE-PUFF 


TIME  1.000/10 
LflXWEN 


riME  1 -QQO'IO 
LflXWEN 


RADIUS  *10  (CM) 
Figure  H-A.  VF-IAX-WENDROFF 


SOLUTION 


SlfiS 


SOLUTION  TINE  l  .000*10 


IS  *10  l CM) 
l-C.  VE-EXACT 


iTITi 


SOLUTION 


RADIUS  >10  ICNJ 
Figure  1 1 -D.  PD-FXACT 


<*31  Ot< 


tlili 


'E-PUI 


iTiTi 


sure  tl-E.  VE-CXACT 


TINE  1.000*10 
PUFF  66 
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3.  TEST  PROBLEM  SCTP-III 
a.  The  Exact  Solution 

In  this  problem,  a  piston  proceeds  with  a  constant  acceleration  into  a 
gas  initially  at  rest  (by  "gas  initially  at  rest"  is  meant  that  the  initial 
conditions  are  as  follows:  velocity  is  zero;  density,  pressure  and  all  other 
fluid  parameters  are  constant).  This  forms  what  is  called  a  compression  wave. 
At  time  tg  =  2Cr/a(y+l),  a  shock  wave  is  formed  (tg  =  time  of  shock  formation, 
Cr  ■  sound  speed  of  the  gas  at  rest,  a  =  acceleration  of  the  piston).  Until 
time  tg,  the  variables  are  continuous  and  the  solution  is  easily  found.*  More¬ 
over,  except  for  one  point  (th«  front  of  the  compression  wave),  the  variables 
are  smooth  prior  to  tg.  The  cotrpression  wave  front  up  to  time  tg  is  Xc(t)  = 
Crt.  After  that  time,  the  compression  wave  front  is  a  shock,  i.e.,  there  is  a 
discontinuity  in  pressure,  density,  velocity,  etc. 

The  solution  for  the  velocity  is 


and 

v(X,t)  -  0  for  X  >  Xc 


*  See  K.  0.  Friedrick's  paper  in  1948  Communications  Pure  and  Applied  Mathematics, 
page  211,  for  an  investigation  of  the  solution  after  shock  formation. 
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Then  the  simple  wav'  formulas  yield 


0  ■  C'(I+^  k) 


p  -  p 


ft) 


2 

Y-l 


P  -  P 


r(‘r) 


2Y 

Y-l 


Notice  that  at  time  tg  «  2Cr/a(Y+l)  and  position  Xg  «  Crtg, 

lim  vx(X»tg)  -  -  ® 

X+X- 


This  indicates  that  a  shock  forms  at  (Xg,tg).  For  further  details  see  Hydrocode 
Test  Problems,  AFWL-TR-67-127. 


The  necessary  data  for  this  problem  are: 

Initial  values:  Pr,  pr,  vr. 

Boundary  values:  At  the  piston  position  Xp  ■  at2  the  velocity  is  vp  *  at. 
There  are  two  variations  of  this  problem: 

SCTP-III-A: 

Pr  «  104  dynes/cm2 
Pr  -  10~6  gm/cm3 
vr  ■  0 

Cr2  «  YPrVr  *  1.4  x  1010  cm2 /sec2 

a  »  Cr/1  sec 

AX  »  10  meters 
Xq  ■  1500  meters 

This  problem  is  run  to  1  second.  The  shock  forms  at  .833...  second. 
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SCTP-III-B 

We  vary  this  from  A  by  setting 

a  =  10  Cr/1  sec 

AX  =  1  meter 
Xq  =  150  meters 

This  problem  is  run  to  .1  second.  The  shock  forms  at  .08333  second. 

b.  The  PUFF  Solution 

The  main  error  made  by  PUFF  is  a  slight  overround  at  X^.  See  Tables  and 
Figures  III. 

c.  The  LAX-WENDROFF  Solution 

The  most  noticeable  error  made  by  the  LAX-WENDROFF  scheme  is  the  oscilla¬ 
tion  just  left  of  X^..  See  Tables  and  Figures  III.  The  time  step  factor  used  was 
.78,  the  artificial  viscosity  factor  used  was  .5. 
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exact  solution  time 


urellt-A.  PD-EXACT 


EXACT  SOLUTION  TIME  8.333‘IG 
SCTP-I I I-fl 


RflOIUS 


TIME  8.333*10 
LflXWEN 


Figure  1 1 1 -A.  VEH.AX-WENDROFF 


8.333*10 


RADIUS  >10  ICR) 
Figure  Ili-B.  PD-EXACT 


SOLUTION 


Figure  I 


0»10 
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SCTP-III-e  PUFF  66 
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Figure  PO-IAX-WENDROFF 
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4.  TEST  PROBLEM  SCTP-IV 
a.  The  Exact  Solution 

In  this  problem  a  piston  has  a  constant  acceleration,  a,  away  from  a  gas 
at  rest.  The  piston  is  traveling  to  the  left.  Eventually,  the  piston  speed 
exceeds  the  speed  with  which  the  gas  can  follow.  This  speed  is  2C  /(y-1),  where 
Cr  is  the  sound  speed  of  the  gas  at  rest.  This  speed  is  called  the  escape  speed 
of  the  gas  and  when  the  piston  exceeds  this  speed  a  vacuum  occurs  between  the 
piston  and  the  edge  of  the  gas.  The  piston  pulls  away  from  the  gas  at  time 

fcv  ~  I a I  (y~1)  and  position  Xp(tv)  »  Xy  =  ^  at2.  solution  for  the  velocity 

is 

v(x>t) . 


for  Xp  <  X  <  Xc,  0  <  t  <  tv,  Xp  -  h  at2,  Xc  -  Crt 
v(X, t)  «  0  for  X>XC  for  all  times. 

For  t  >  ty  there  is  a  vacuum  from  Xp  to  the  gas  front.  The  gas  front  is  at 
2cr 

\  ~  Therefore  the  velocity  is  really  meaningless  in  this  region. 

But  the  pressure,  density,  and  sound  speed  are  all  zero  in  this  region.  From  the 
gas  front  position  on  to  X^  the  above  formula  for  the  velocity  holds.  Once  the 
velocity  is  known  the  simple  wave  formulas  yield 


C 


P 


P  »  P 


r 
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The  necessary  data  for  this  problem  are: 

Initial  values:  Pr,  pr,  vr 

Boundary  values:  At  the  piston  position  Xp  =  ijat2  the  velocity  is  vp  --  at 
There  are  two  variations  of  this  problem: 

SCTP-IV-A: 

Pr  =  10**  dynes /cm2 
Pr  *  10-6  gm/cm3 
vr  =  0 

C2  =  yPrVr  =  1.4  x  10 10  cm2 /sec2 
a  =  -  Cr/1  sec 
AX  =  1 )  meters 
Xq  *  1500  meters 

This  problem  is  run  to  10  seconds;  the  vacuum  forms  at  5  seconds. 

SCTP-IV-B: 

Same  as  A  but  a  «  -  10Cr/l  sec,  run  the  problem  to  1  second  and  the  vacuum  occurs 
at  .5  second. 

b.  The  PUFF  solution 

PUFF's  main  errors  in  this  problem  are  at  X^  and  Xp.  At  Xc  the  error 
PUFF  makes  is  an  underround  in  pressure,  density,  velocity,  and  internal  energy. 
The  error  PUFF  makes  just  to  the  right  of  Xp  is  due  to  the  fact  that  PUFF  is  a 
Lagrangian  code  and  the  mass  that  was  originally  in  a  zone  remains  in  the  zone 
and  therefore  the  density  and  pressure  can  never  go  to  zero.  See  Tables  and 
Figures  IV. 

c.  The  LAX-WENDROFF  Solution 

Since  a  vacuum  occurs  in  this  problem  and  the  LAX-WENDROFF  code  uses 
specific  volume  as  a  variable  it  cannot  run  this  problem. 


ERRORS  ON  SCTP-IV-A 


Int.  Energy _ Sum  Kin.  Energy  Sum  Tot.  Energy 
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5.  TEST  PROBLEM  SCTP-V 
a.  The  Exact  Solution 

This  is  called  the  shock  tube  problem.  It  is  an  example  of  the  more 
general  Riemann  problem.  The  Riemann  problem  is  that  of  determining  the  flow 
after  the  conjunction  of  two  states,  left  state  and  right  state,  with  P?>  p^,  v^ 
the  constant  values  of  the  left  state  and  Pr,  pr,  vr  the  constant  values  of  the 
right  state.  In  the  shock  tube  problem,  vr  and  Vj  are  no  longer  arbitrary  but 
are  set  to  zero.  So  the  problem  may  be  interpreted  as  the  determination  of  the 
flow  after  removal  of  the  membrane  separating  two  constant  states  at  rest.  As 
a  convention,  take  P^  >  Pr  £  0.  Then  in  the  code  test  problem  the  three  poss 
bilities,  >  pr,  p^  “  ox,  and  p^  <  pr  will  be  explored. 

At  time  zero,  the  membrane  is  removed.  The  resultant  action  is  a  rare¬ 
faction  wave  traveling  into  the  left  3tate  and  a  shock  traveling  into  the  right 
state.  The  velocity  is  v^  =  0  to  the  left  of  the  rarefaction  wave.  From  the 
left  of  the  rarefaction  wave  to  the  right,  the  velocity  rises  linearly  from  0  to 
vm  >  0.  The  velocity  is  constant  at  vffl  from  the  right  of  the  rarefaction  wave 
rightward  toward  the  shock.  Ac  the  shock,  the  velocity  jumps  from  vn  >  0  down 
to  vr  ■  0.  The  pressure  drops  continuously  across  the  rarefaction  wave  from  P^ 
to  Pm.  The  pressure  has  the  value  Pm  constantly  from  the  right  of  the  rarefaction 
wave  to  the  shock.  The  pressure  drops  from  Pm  to  Pr  across  the  shock.  The  densi¬ 
ty  drops  continuously  across  the  rarefaction  wave  from  p^  to  a  value  p  which 
it  maintains  from  the  right  of  the  rarefaction  wave  to  the  point  in  the  fluid 
where  the  initial  discontinuity  was  and  there  the  density  jumps  up  to  Pmr. 

From  the  initial  discontinuity  point  to  the  shock,  the  density  jumps  down 
from  cmr  to  pr,  which  value  is  maintained  all  the  way  to  the  right. 

The  left  side  of  the  rarefaction  wave  is  at  Xc(t)  *  X<,(0)  -  C^t,  where 
Xs(0)  is  the  position  of  the  shock  at  time  zero  which  is  also  the  position  of  the 
initial  discontinuity. 

The  right  side  of  the  shock  is  at 

XR(c>  -  XS(0)  -  (c  -  lii  vm)t 
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The  shock  wave  is  at 


Xs(t)  =  X3(0)  +  vst 
where  Vg  is  the  velocity  of  the  shock 


vs  -  vr  +  ‘■r) 

The  initial  discontinuity  point  of  the  fluid  is  at  X0(t)  =  Xg(0)  +  v| 
middle  values  vm  and  Pm  are  determined  by  simultaneously  solving 


and 


vm  =  vr  +  (Pm  -  P,.)  J 
m  r  '  m  r7  ^ 


2V. 


(V+I)  Pm  +  (y-1)  Pr 


t.  The 
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Solution  Summary: 


LEFT 

REGION 


For  X  <  X^(t),  the  values  are  P.,  p  ,  v 

X  X  X 


For  Xc(t)  <  X  <  XR(t),  the  values  are 


v(X,t) 


X  -  xc(t) 

XR<fc)  ~  XC(t)  Vm 


RAREFACTION 

REGION 


C  -  C, 


(‘-■n) 


P  »  p 


2Y 

-ft)3 


P  *  P 


ft)” 


For  XR(t)  i  X  <  Xc(t),  the  values  are  P  ,  v 
o  in  ui 


MIDDLE 

REGION 


For  XR(t)  <  X  <  XD(t),  the  density  is  pm^ 


For  XD(t)  <  X  <  Xg(t),  the  density  is  p 


RIGHT 

REGION 


For  X  >  Xg(t),  the  values  are  Pr,  pr,  v„ 
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The  necessary  data  for  this  problem  are 

INITIAL  VALUES:  Pr,  pr,  vr,  P£,  p£>  v£ 

BOUNDARY  VALUES:  At  X  =  0,  hold  the  values  at  P£,  o->  and  Vj,  at  Xq 
(the  right  boundary)  hold  the  values  at  Pr,  ,  vr. 

There  are  three  variations  of  this  problem 

SCTP-V-A: 

Xg(0)  =  100  meters 
AX  =  1  meter 
P£  =  108  dynes /cm2 
p£  =  10-5  gm/cm3 
v£  =  0 

Pr  =  104  dynes/cm2 

Pr  *  10-6  gm/cra3 

v  -  0 
r 

Xq  *  250  meters 

These  values  imply  the  following  values: 

Pm  =  1.888  x  107  dynes/cm2 
vm  =  3.964  x  106  cm/sec 
om£  =  3.040  x  10~6  gm/cm3 
Pmr  =  5.982  x  10"6  gra/cm3 
Vg  i  4.760  x  108  cm/sec 
This  problem  was  run  to  2  x  10" 3  seconds. 

SCTP-V-B: 

This  problem  is  the  same  as  SCTP-V-A  except 
Xs(0)  =  250  meters 
p£  *  10" 6  gm/cm3 
Xq  *  500  meters 
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These  values  imply  the  following  values: 

Pm  =  4.610  x  10  7  dynes/cm2 
vm  *  6.196  x  106  cm/sec 
Pmi  *  5.751  x  10“  7  gm/cm3 
omr  ^  5.992  x  10"6  gm/cm3 
vs  i  7.437  x  106  cm/sec 

SCTP-V-C: 

This  problem  is  the  same  as  A  except 
Xg(O)  *  250  meters 
*  10“ 6  gm/cm3 
Pr  -  10“5  gm/cm3 
Xq  *  500  meters 

These  values  imply  the  following  values: 

Pm  *  7.406  x  10 7  dynes /cm2 
vm  i  2.484  x  106  cm/sec 

Pm£  A  8.070  x  10~7  gm/cm3 

pmr  i  5’995  x  10-5  g“/cm3 
Vg  *  2.981  x  10  6  cm/sec 

b.  The  PUFF  Solution 

On  SCTP-V-A,  the  most  noticeable  error  was  a  smearing  of  the  density 
discontinuity  at  Xrj.  The  only  errors  were  the  typical  underrounds  at.  i  overrounds 
at  corners. 

On  SCTP-V-B,  there  was  a  bit  of  oscillation  in  the  density  in  the 
compressed  region  and  a  little  overshoot  in  velocity  and  an  undershoot  in  pressure 
at  XR. 

On  SCTP-V-C,  the  dominant  error  was  a  siight  undershoot  in  the  pressure 
at  XR.  The  other  error  was  a  slight  undershoot  in  the  pressure  at  XR. 

For  more  details  see  Tables  and  Figures  V. 
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c.  The  LAX-WENDROFF  Solution 

In  order  to  run  this  problem  it  was  found  necessary  to  cut  the  first  time 
step  and  artificial  viscosity  factor  down  to  one-twentieth  the  normal  time  step, 
cut  the  second  time  step  and  artificial  viscosity  down  to  two- twentieths  of  the 
normal  time  step,  etc.,  until  the  twentieth  step  and  thereafter  allow  the 
normal  time  step  and  artificial  viscosity  factor.  The  time  factor  used  was  .78 
and  the  artificial  viscosity  factor  used  was  .5. 

Ihe  most  noticeable  difference  between  PUFF  and  LAX-WENDROFF  in  SCTP-V-A 
is  the  pronounced  spikes  at  X,,  and  Xg  in  LAX-WENDROFF  (see  Figures  V-A) . 

In  SCTP-V-B  the  spikes  are  nor  so  bad  bu;  there  is  quite  an  oscillation 
in  the  velocity  just  right  of  X^  (see  Figures  V-B). 

In  SCTP-V-C  the  overshoot  in  the  velocity  at  X^  has  grown  mere  pronounced 
(see  Figures  V-C). 
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Sum  Int.  Energy  Sum  Kin.  Energy  Sum  Tot.  Energy 
EXACT  2.17716  x  1012  5723215  x  10 11  2.50037  x  1012 
KjFF  2.17807  x  1012  3.21760  x  10^  2.49983  x  1012 


LAX-WENDROFF 


2 


3 


EXACT  SOLUTION  T I  HE  2.QQ0>10 


KflOiUS  >10  (CM) 


Tine  2 .GOQ'IQ 

PUFF  66 


RflOI'JS  *10  (CM) 


:  363  TIME  2.QD0>10 

-V-fl  LflXWEN 


N3MXB1 
0lv0Q0*2  3WU 


Figure  V-A.  VE-LAX-WENDROFF 


[XflllT  SOLUTION  TIME  2-000»10 
SCTP-V-8 


PO-EXACT 


EXACT  SQLUT I  UN  TIME  2.00Q>ICI 


RADIUS  '10  ( CM l 


527 


1527 


EXACT  SOLUTION  TIME  2- OOO’IO 
5CTP-V-C _ 


b 

% 

i 


Figure  V-C.  VE-EXACT 


r  I  ME  2  -QOO’IO 
PUFF  66 


PD-PUFF 


iliTI 


IOBH 


“YCl.t  337  T  l  HE  2.000M3 

SCTF-v-f.  LFXWEN 
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6.  TEST  PROBLEM  SCTP-VI 
a.  The  Exact  Solution 

This  problem  is  the  collision  of  two  shock  waves.  It  is  another  special 
case  of  the  Riemann  problem.  Proceeding  from  left  to  right,  the  initial  values 
are  P^»  p^,  v^  connected  by  a  right-facing  shock  to  PQ,  pf, ,  v^,  which  in  turn  is 
connected  by  a  left  facing  shock  to  Pr,  pr,  vr.  As  a  convenient  convention  one 
takes  >  Pr  1  Pq. 

After  collision  a  shock  travels  back  to  the  la£Z  &aP  '  3hock  travels  on 
to  the  right  from  a  middle  region  in  which  the  velocity  and  ueessure  are  constants 

vm  and  V 

For  a  shock  facing  to  the  right 


v«  "  vr  +  *r(r.) 

and  lot  a  shock  facing  to  the  left 

■  vi  -  «*(*■) 

where 

2V 

da(p)  -  (p  -  Pa)  (V+1)?  +a(-r-i)Pa 

From  these  two  equations,  Pffi  and  vm  are  determined. 

The  density  profile  proceeding  from  left  to  right  is  p^,  then  it  jumps  up 

to  om«  at  the  left-facing  shock,  then  at  the  point  of  collision  (in  the  Lagracgian 

coordinates)  the  density  jumps  to  pmri  then  at  the  right-facing  shock  the  density 

jumps  down  to  pr.  The  Rankine-Hugoniot  relation  determines  p  s  and  p  .  That  is, 

©~  tnr 


0  •  w 


Pj  +  Pr 
2 
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+  P. 


e”r ' er + (v  -  vr) 


where 


*nd,  of  course, 


?«*  •  p;J.  vmr  .  p;i 


of  the  shor'c  veloM  rim*,  —  _  , 

■tocl ties  ney  be  computed  by 


VS  "  v£  +  vr  +  v£  +  v  \2  +  c  2 

A  r/  r 

!"“* +  “  “k,n  £or  ,hock,  Md .  la 

For  ersnple,  prior  to  collision  (let  t  .  *ft~fsclng  shocks. 

04  £h*  **«*«“*  -*  the  the  Veloclty 

v?i*^v‘+i?v0+ 


*°d  ”*l0CUy  °f  *•  WMd,  shock 


the  shock  on  the  right)  is 


vSr.i|*Vo+¥vr_^Vo+¥vJ  +  Cr2 

2:  rrn  “r >  w  the  °f  -  — .  <,... 


V  *  1+Y  3-y 

VSr  m  —  v£  +  ^  Vjt  - 


(it!  v  +  3^r  \2 

\a  vi+-4~v»;  +v 
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and  the  velocity  of  the  right-facing  shock  (i.e.,  the  shock  on  the  right)  is 


where  as  before  C  stands  for  the  isentropic  sound  speed. 
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Solution  Summary: 

Prior  to  collision  (t  <  tco£) 


LEFT 

REGION 


For  X  <  Xg£(t)  -  Xg^CO)  +  vg£t,  the  values  are  P£,  p£,  v£ 


MIDDLE 

REGION 


RIGHT 

REGION 


j  For  XSi^  <  x  <  xsr(t)  *  xsr(°)  +  vsrt,  the  values  are  PQ, 
|  For  X  >  X8r(t),  the  values  are  Pr>  pr,  vr 


After  collision  (t  >  tco£) 

LEFT  f  *  * 

REGION  |  For  X  <  X^(t)  -  Xco£  +  vg£(t-tco£),  the  values  are  P£,  p£,  v 


(  For  4<«  <  X  <  XSr(t>  ■  Xcot  +  W. 


MIDDLE 

REGION 


the  values  ate  P„,  v,  and  P|jl  for  X  <  Xcol  +  v„(t-tC0J[) 
and  omr  for  X  >  X.0.  +  vm(t-t  .) 


RIGHT 

REGION 


For  X  >  Xsr(t),  the  values  are  Pr>  pf,  vf 
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The  necessary  data  for  this  problem  are: 

INITIAL  VALUES:  P£,  PQ,  p0,  vQ,  Pr 

From  these  all  other  initial  values  are  determined. 

BOUNDARY  VALUES:  At  X  *  0  (the  left  boundary)  hold  the  values  at  Pft,  P£, 
V£.  At  Xq  (the  right  boundary)  hold  the  values  at  Pr, 
Pr.  vr. 

There  are  two  variations  of  this  problem: 

SCTP-VI-A: 

AX  *  1  meter 
XS£(0)  =  75  meters 
XSr(0)  =  125  meters 
Xq(0)  =  200  meters 
PQ  **  104  dynes /cm2 
pQ  *  10-6  gm/cm3 
P^  *  10 8  dynes /cm2 
Pr  =  10 7  dynes /cm2 
vQ  *  0 

These  values  then  determine  the  following  values: 

p^  *  5.997  x  10  6  gm/cm3 

Pr  "  5.97  x  10-6  gm/cm3 

v^  A  9.13  x  106  cm/sec 

vr  *  -2.88  x  106  cm/sec 

v^  *  1.095  x  107  cm/sec 

vSr  A  -3.46  x  106  cm/sec 

tC0£  *  3.468  x  10-4  sec 

X  A  1.13  x  104  cm 
col 

Pm  A  3.66  x  108  dynes/cm2 
Pm£  A  i,43  x  IQ-5  gm/cm3 
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Pmr  ■  3.09  x  10“ 5  gm/cm3 
vm  i  1.96  x  106  cm/sec 
VS£  *  3.75  x  105  cm/sec 
vgr  i  5.72  x  106  cm/sec 
This  problem  was  run  to  7  x  10"4  sec. 

SCTP-VI-B: 

This  problem  is  the  same  as  A,  except  ■  Pr  ■  108  dynes/cm2. 

This  yields  the  following  values: 

p£  *  pr  *  5.997  x  10“6  gm/cm3 
«  -vr  A  9.13  x  106  cm/sec 
vs  »  “vsr  *  1*095  x  107  cm/sec 
tcoS.  “  2.282  x  10“4  sec 

Xco£  -  1.00  x  104  cm 
Pn  *  7.995  x  108  dynes /cm2 

vm  -  0 

pm£  *  pmr  i  2*098  x  10~S  go/cm3 

*  *  .  _  , 

“VS£  *  vSr  “  3.65  x  10b  cm/sec 

This  problem  was  run  to  7  x  IQ”4  sec. 

b.  The  PUFF  Solution 

The  major  errors  in  evidence  were  the  spikes  in  the  density  and  internal 
energy.  Hot-thin  spikes  resulted  from  the  initial  discontinuities  and  a  cold- 
thick  spike  resulted  from  the  shock  collision.  For  more  details  see  Tables  and 
Figures  VI. 

c.  The  LAX-WENDROFF  Solution 

In  addition  to  the  spikes  observed  in  the  PUFF  solution  there  is  also,  a 
considerable  amount  of  oscillation  in  evidence  in  the  LAX-WENDROFF  Solution.  The 
time  factor  used  was  .39  and  the  artificial  viscosity  factor  used  was  .25;  as  in 
SCTP-V  the  time  and  viscosity  factors  were  multipled  by  one-twentieth  on  the  first 
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fHrr'v^xtFr •~1\*J > r?<xr ~ rv?4 6 . v "v~r * ^zn-i ; iT&ttT^T' v<-"<iJ 


time  step,  two- twentieths  on  the  second,  etc.,  until  the  twentieth  time  step  and 
thereafter  when  they  were  left  at  the  values  of  .39  and  .25  respectively.  For 
more  de tells  see  T  .b'  "'j  ae,J  ~l  cres  VI. 
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7.  TEST  PROBLEM  SCTP-VII 
a.  The  Exact  Solution 

In  this  problem  one  shock  wave  overtakes  another.  This  problem  is 
another  special  case  of  the  Riemann  problem.  Two  shock  waves  are  traveling  in 
the  same  direction,  which  is  taken  to  the  right.  When  two  shock  waves  are 
traveling  in  the  same  direction,  the  one  behind  will  always  overtake  the  one  in 
front.  After  overtake  time,  a  rarefaction  travels  back  to  the  left  (for  y  <  5/3) 
and  a  stronger  shock  travels  on  to  the  right  and  there  is  a  middle  region  in 
which  the  velocity  and  pressure  are  constants  vm  and  Pm. 

Proceeding  from  left  to  right,  the  initial  values  are  P^,  v^  connected 
by  a  right-facing  shock  to  P^r,  pfj.  which  in  turn  is  connected  by  a  right-facing 
shock  to  P_,  pr,  vr. 

After  overtake  vm  *  vr+<j>r(Pni)  for  the  shock  traveling  to  the  right  and 
vm  *  v£  -  ^(Pn,)  for  the  rarefaction  traveling  to  the  left.  Recall  that 


*r 


3!e 

+ 


(Y-l)P 


r 


and 


From  the  above  relations  v  and  P_  are  determined. 

m  m 

In  the  middle  region  there  will  be  two  values  for  the  density;  P£m  - 
to  the  left  of  the  overtake  point  and  pmr  *  V”£  to  the  right  of  the  overtake 
point.  The  Rankine-Hugoniot  relation  determines  Vmr  by 


0 


mr 


•r  »,) 


p»  +*r 
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and  is  determined  from  the  fact  that  the  entropy  does  not  change  through  a 
rarefaction;  therefore, 


Let  vg£,  Xgi  and  Vgr,  Xgr  be  the  velocities  and  positions  of  the  left  and 
right  shocks  prior  to  overtake;  (XQ,  tQ)  be  the  point  where  overtake  occurs;  Vg, 

Xc  be  the  velocity  and  position  of  the  shock  after  overtake;  Xr  be  the  left  side 
of  the  rarefaction  wave;  Xr  be  the  right  side  of  the  rarefaction  wave  and  XD  be 
the  position  of  the  point  in  the  fluid  where  overtake  occurs. 

( 


f 


i 
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Solution  Summary: 

Prior  to  overtake  (t  <  tn) 


LEFT  ) 

REGION  | 

MIDDLE  ( 

REGION  | 


For  X  <  Xg£(t),  the  values  are  P£,  v£,  p£ 


RIGHT 

REGION 


\  For  XS£(t)  <  x  <  xsr(t)*  the  values  are  P£r,  v£r>  P£r 
|  For  X  >  Xgr(t),  the  values  are  Pr,  vr,  pr 


After  overtake  (t  >  tQ) 


REGION 


MIDDLE 

REGION 


For  X  <  X0  +  (v£~c£) (t-tQ)  -  Xc(t),  the  values  are  P.,  v.. 


£>  V£,  p^ 


For  xc(r)  <  x  <  x0  +  (-  c,  +  0±i  v„  +  iii  v,)  <t-t0,  .  xR(t>, 
the  velocity  goes  linearly  from  v£  at  Xc(t)  up  to  vm  at 
XR(t)  and 


C  -  0£  -  (v£-v) 


P  -  P, 


(*r 


P£  "  P* 


(0 


For  the  region  Xs(t>  <  X  <  X„  +  v„(t-t„>  -  XD!t),  the  values  are 

Pm»  0£o’  and  vn.  For  the  .region  Xpit)  <  X  <  X0  +  vs(t-t0)  «  Xs(t), 
the  values  are  P^,  par,  va. 


RIGHT 

REGION 


For  X  >  Xg(t),  the  values  eze  Pr,  pr,  vr. 
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The  necessary  data  for  this  problem  are: 

INITIAL  VALUES:  P£,  ?lr,  ?r,  c*r, 

From  these  values  all  other  initial  values  are  deternined. 


BOUNDARY  VALUES: 


At  X  *  0  (the  left  boundary),  hold  the  values  at  P£, 
pj.,  v^  and  at  X-  (the  right  boundary)  hold  the  values 


at  Pr,  er. 


Now  the  specific  numerical  values  are  presented  for  SCTP-VII: 


AX  -  1  meter 
Xst(0)  *  25  aeterr- 

XSr(0)  *  100  meters 

Pr  *  104  dynes /cm2 

Pr  «  10"6  gn/ca3 

vr  *  0 

P£r  »  108  dynes/cm2 
P^  *  1012  dynes/cxsz 
Xq  *  200  meters 
These  values  yield 

Ci  *  1.97  x  108  co/sec 
v.  A  3.82  x  10s  ca/sec 
p£  A  3.60  x  10~5  g«/cm3 

vsi  *  4*56  x  i®8  c*/sec 

vSr  A  1.095  x  107  cn/sec 
tQ  A  1.683  x  10“  5  sec 
X0  A  1.018  x  104  cm 
Pu  A  3.32  x  1011  dynes/ca2 
v£r  A  9.13  x  108  cm/sec 
PU  A  5.997  x  10'6  ga/ca3 
vQ  A  5.26  x  106  ca/sec 
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Vg  A  6.31  108  cm/sec 

pmi  1.63  x  10" 5  gto/cm3 

pmr  “  6.000  x  10“6  gm/cm3 

This  problem  was  run  to  3  x  10~5  sec. 

b.  The  FUFF  Solution 

As  in  SCTP-VI,  the  major  errors  in  evidence  were  the  spikes  in  density 
and  internal  energy.  Hot-thin  spikes  resulted  from  the  initial  shock  disconti¬ 
nuities  and  a  cold-thick  spike  resulted  from  the  shock  overtake.  For  more  details, 
see  Table  and  Figures  VII. 

c.  The  LAX-KENDROFF  Solution 

In  addition  to  the  spikes  observed  in  the  PUFF  solution,  there  is  also 
more  oscillation  in  the  LAX-WENDROFF  solution.  The  time  factor  used  was  .39,  the 
artificial  viscosity  factor  used  was  .25,  and  both  factors  were  multiplied  by 
one-twentieth  on  the  first  time  step,  two-twentieths  on  the  second,  etc.,  until 
the  twentieth  time  step  and  thereafter  when  they  were  left  at  the  values  of  .39 
and  .25  For  mora  details,  see  Table  and  Figures  VII. 
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SECTION  III 
CONCLUSIONS 


The  most  apparent  difference  between  the  PUFF  and  LAX-WENDROFF  solutions  is 
the  greater  tendency  of  the  LAX-WENDROFF  scheme  to  oscillate. 

In  those  flows  in  which  there  are  no  strong  shocks  of  strong  rarefactions  or 
vacuums,  the  LAX-WENDROFF  scheme  is  more  accurate  than  PUFF.  However,  in  those 
flows  in  which  there  are  strong  shocks  or  strong  rarefactions  or  vacuums  the  PUFF 
scheme  is  more  accurate.  The  LAX-WENDROFF  scheme  cannot  handle  vacuums  because 
of  the  use  of  the  specific  volume  instead  of  the  density  as  a  fluid  variable.  It 
appears  that  the  LAX-WENDROFF  scheme  could  be  improved  by  using  an  artificial 
viscosity  of  the  type  used  in  PUFF.  And  in  general  it  appears  that  it  might  be 
possible  to  combine  the  better  features  of  PUFF  and  LAX-WENDROFF  to  produce  a 
superior  hydrocode.  This  will  be  investigated  and  discussed  in  a  forthcoming 
report. 
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